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Abstract 

We propose an evolutional scenario of the universe which starts from 
quantum states with conformal invariance, passing through the inflation- 
ary era, and then makes transition to the conventional Einstein space- 
time. The space-time dynamics is derived from the renormalizable higher- 
derivative quantum gravity on the basis of a conformal gravity in four 
dimensions. Based on the linear perturbation theory in the inflationary 
background, we simulate evolutions of gravitational scalar, vector and 
tensor modes, and evaluate the spectra at the transition point located at 
the beginning of the big bang. The obtained spectra cover the range of 
the primordial spectra for explaining the anisotropies in the homogeneous 
CMB. 
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1 Introduction 



After passing many theoretical as well as experimental tests, the general 
theory of relativity has been established as the fundamental theory of gravity 
capable of describing the universe. On the other hand, tracing the history of the 
universe, the space-time in its early epoch would be totally fluctuating quantum 
mechanically so that geometry lose its classical meaning. This would imply that 
there is a transition from quantum space-time to classical space-time, and there 
should exist a dynamical scale separating these two phases. 

There is a possibility to observe the instance of the transition, because we 
can trace the past guided by known physical laws as far as the classical general 
relativity holds. Valuable information on physical processes taking place in the 
expanding universe has been recorded in the cosmic microwave background ra- 
diation (CMB) as tiny anisotropics. Angular power spectra of the anisotropies, 
recently observed by the Cosmic Background Explorer (COBE) [T] and the 
Wilkinson Microwave Anisotropy Probe (WMAP) [51 [3], are, roughly speak- 
ing, projection of the history of universe for the period between its birth to 
the present. It gives us an amazing hope that if we believe the idea of infla- 
tionary universe [4} meaning an extremely rapid expansion without global 
thermalization, the long-distance correlation in observed anisotropies can pro- 
vide information about dynamics of the period before the universe grew to the 
Planck scale. We are now at the stage of revealing and verifying the quantum 
aspect of universe. 

The model of space-time transition proposed in this paper emerges from the 
renormalizable higher-derivative quantum theory of gravity developed on the 
basis of the conformal gravity in four dimensions. In this theory, at very high 
energies beyond the Planck scale, quantum fluctuations of the conformal mode 
in the metric field are dominated, and it is treated non-perturbatively. The 
space-time is described by a conformal field theory whose dynamics is governed 
by conformal invariant gravitational actions. The conformal field theory loses 
its validity at the dynamical scale indicated by the asymptotic freedom of the 
unique dimcnsionless coupling constant introduced for the traceless tensor mode. 
At about this point the universe is expected to make a transition from the 
quantum space-time to the classical Einstein space-time. 

There are three mass scales in the model, namely the Planck mass Mp, 
the dynamical scale Aqq, and the cosmological constant Acos- We set their 
ordering as 

Mp» AQG»A^(^g. (1.1) 

We shall obtain an evolutional scenario according to the order starting infla- 
tion driven by quantum effects of gravity without adding any artificial field by 
hands [5]. The inflationary model induced by quantum effects of gravity was 
first proposed by Starobinsky _5j . We here develope along the idea and propose 
an evolutional scenario of the universe summarized as follows: the conformal 
symmetry begins to be broken about the Planck scale to form an inflationary 
universe with the expansion time constant of order of the Planck mass, and 
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completely broken at the dynamical scale turning to the classical Einstein uni- 
verse. Then, energies stored in extra degrees of freedom in higher-derivative 
gravitational fields shift to matter degrees of freedom, causing the big bang. 
The primordial power spectra obtained from the two-point correlation func- 
tions of gravitational fields show a significant character of the transition which 
is expected to be observed cosmologically. 

The aim of this paper is to clarify why we can observe the Planck scale 
phenomena today from the tiny CMB anisotropics. Evolution of the universe 
is described by the equations of motion taking effects of running coupling as a 
time-dependent function. We find an inflationary solution that starts from the 
Planck scale and end at the dynamical scale where the transition of space-time 
occurs. Since the inflationary homogeneous solution is stable, the fluctuations 
from this solution will be slowly diminishing during the inflationary era, and thus 
the linear perturbation about the inflationary solution becomes applicable. We 
evaluate gravitational fluctuations perturbed about the homogeneous solution: 
the scalar perturbation so-called Bardeen potential, the vector perturbation 
and the tensor perturbation. We obtain the spectra at the transition point, 
which should be used for the primordial spectra to analize the data observed by 
WMAP. The rest of the paper organized as follows: we summarize the model 
of quantum gravity in section 2 and construct an evolutional scienario of the 
universe in section 3. In section 4 we study the linear perturbation theory in 
the inflationary background. We then simulate evolutions of perturbations and 
obtain primordial spectra in section 5. We conclude in section 6. 

2 Renormalizable Quantum Gravity 

Since it has been recognized that any attempt to quantize Einstein grav- 
ity perturbatively [H [H [9] cannot be succeeded [TOj, most of researchers in 
this field feel the necessity of quantizing gravity in either a non-perturbative 
way or others. Historically, the renormalization problem is tackled introduc- 
ing four-derivative terms to the Einstein-Hilbert action [TTl [121 US] , so that the 
gravitational coupling constant becomes dimensionless, and at the same time 
we can avoid the unboundness problem. Among various models, we employ the 
conformal gravity without term [Ml [15] . 

The Model The conformal gravity is defned by the action, 

/( 1 A'P 
d*xy^}^--Cl,^, - + -fR - Acos 

-\Tr - i i^d^Xd.X + IrX^^ + . . j, (2.1) 

where we write the reduced Planck mass and the cosmological constant as 
Mp = I/vottG and Acos, respectively. There are two conformal invariant 
four derivative actions in four dimensions: the square of the Weyl tensor C^i,^^ 
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and the Euler term G4. The Weyl tensor is the field strength for the traceless 
tensor mode, and t is the dimensionless couphng constant. The constant b is 
introduced to renormahze divergences proportional to the Euler term, which is 
not an independent constant because it does not have a kinetic term. 

The third term is the Einstein-Hilbert action, which will dominate at energies 
below the Planck scale. Ff^^, is a gauge field strength, and X is a scalar field 
with conformal coupling. In addition to these fields, we may consider fermions, 
which are conformally invariant in any dimensions, without mass terms. 

The beta function of the renormalized coupling constant tr indicates asymp- 
totic freedom. This implies that at high energies beyond the Planck scale, the 
coupling constant gets small, and accordingly configurations with the vanishing 
Weyl tensor will dominate. Then, the singular configulation such as a black 
hole at which the Riemann-Christoffel curvature tensor is divergent is excluded 
quantum mechanically. This is a prominent feature of the renormalizable quan- 
tum gravity based on the conformal gravity, in contrast to the theory based 
on the Einstein-Hilbert action because such a configuration has the vanishing 
scalar curvature so that its quantum weight in the path integral is unity. 

The partition function is defined by the functional integral over the metric 
field Qf^^. We here decompose the metric field into the conformal mode <j), the 
traceless mode h\ and the background metric 

5m- = (2.2) 

and 

.9m- = = 9,x {6i + + ■■■), (2.3) 

with tr{h) = h\ — 0. Since we will consider a case with large running coupling 
constant, we include the coupling t in the traceless mode h\ instead of th\. 
The signature of the metric is taken as (—1,1,1,1). The contraction of the 
indices of /i^^ is done by using the background metric g^^i,. In the following, 
quantities with the hat and the bar on them are defined in terms of the metric 
(jfj^i, and 5^1^, and the contraction of them is done by using the metric and 
g^,y, respectively. Unless it is not specified, the flat background — rj^^ is 
used, which defines the comoving frame with the coordinate = {"q^x^) and 
= (<9,,,9i)- 

Techniques to treat diffeomorphism invariance are developed in the two di- 
mensional theory at the end of 1980's [161 El IHl HSIj and extended later in 
four dimensions [20l[2ll[22l[23l[2l[Ill[2Sl[26l[l5]. We change the path-integral 
measures from the diffeomorphism invariant measures to the practical mea- 
sures defined on the background metric g^v. In order to recover the diffeomor- 
phism invariance the Wess-Zumino term S |27j related to conformal anomalies 
PSI Uni 1301 131] in the action is necessary as the Jacobian, and the partition 
function is expressed as 

Z - J ^^^^^^ exp {.^(^, ,) + UiA, X,g)}. (2.4) 
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The induced action S contains the kinetic term for the conformal field. At the 
lowest order in the coupling, it is given by the Riegert action [20] 

(2.5) 

where (7A4 is a conformal invariant fourth order self-adjoint operator for a 
scalar field variable and = V^Va. The coefficient 61 has been computed 
within the lowest order as [22 

^r = ^(iVx + ^A^w + 62iV.)+Ig, (2.6) 

where Nx, iV-yv and TVa are the numbers of scalar fields, Weyl fermions and 
gauge fields added to the action, respectively. The dots indicate new vertices 
like (/>""'"^A40, 0"C'^,yA(Ti ^^'^ '/'"^''(-^^ly) induced at higher order of the coupling 

MM- 

Since the conformal variation of the metric g^i, is equivalent to the variation 
of the conformal field, the trace of derived stress tensor must be equal to the 
equation of motion of the confromal field. The second and third terms in the 
Riegert action (|2.5p ensure such a diffeomorphism invariance condition within 
the linear approximation discussed in Section 4. 



Running coupling constant The beta function of the renormalized coupling 
tr for the traceless mode has been calculated to be /3 = —(3oi^ with /3o = {(iVx + 
+ 12iVA)/240 + 197/60}/(47r)2 within the lowest order of perturbation 
[131 fT4] . It indicates the dynamics of the traceless tensor mode asymptotically 
free. The effective action in the momentum space is expanded in tr as 
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CU. (2-7) 



for the Weyl tensor term, where fc is a comoving momentum defined on the flat 
background. The running coupling constant is approximately written as 

/3o log(p2/A^(3) 

where p is a physical momentum defined by p — k/e'^. The coupling constant 
is a measure of the degree of deviation from conformal field theory. The scale 
parameter Aqq — /iexp{— 1/2/30^^} with /i being a renormalization mass scale 
represents the energy scale where the metric field changes its appearance from 
the conformal mode and the traceless mode to the conventional Einstein gravity. 

As discussed later, in the infiationary era the conformal field, </>, grows large 
and the physical momentum diminishes exponentially. Accordingly, the running 
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coupling constant tr{p) gets large during the inflation. As energies become 
lower to the dynamical scale Aqg, the running coupling constant diverges and 
the Weyl terms in the effective action will disappear. There, the conformal 
invariance terminates and the background metric appears as a physical measure 
of space and time. Beyond the energy smaller than this point the Einstein 
gravity dominates the dynamics. 

3 Evolutional Scenario of The Universe 

The renormalizable conformal quantum gravity suggests that there are four 
stages in the evolution of the universe divided by three mass scales. The first 
is quantum phase dominated by the dynamics of the conformal gravity action, 
which gradually turns to the second stage, i.e. the inflationary era. The third 
stage is the Einstein universe mostly ruled by the classical theory of general 
relativity. The last is the present deSitter evolution where the universe expanded 
sufficiently and the cosmological constant becomes effective. The aim of this 
section is to construct an evolutional model connecting the quantum gravity 
phase and the Einstein space-time. 

Inflationary phase In the very early epoch when the space extends much 
smaller than the Planck scale, typical energy is higher than Mp. Corrections 
of order t'^ in the asymptotic free regime can be neglected and the dynamics is 
governed by conformal field theory of gravity. As the universe expands typical 
energy gets lowered to Mp, the Einstein action gradually becomes effective, 
so that the conformal symmetry starts to be broken to develop the classical 
solution. The effective action valid for such an energy regime is given by the 
sum of the Riegert action (|2.5p . the Einstein action and the matter action. 
The homogeneous equations of motion has a stable inflationary solution [Sj . At 
the Aqq scale, this symmetry is completely broken and the universe make a 
transition to the classical space-time. In this section, we consider a dynamical 
effect arising from the effective action under the inflationary background. 

The equations of motion simplified in the local form are obtaind in the way 
that the coefficient of the Riegert action of order t'^ corrections is written as [15] 

bi ^ &i (1 - aitl + •••)= biBoiU), (3.1) 

where ai is a positive constant, and Bq is assumed to be 

Here, k is a parameter for taking into account of the higher order perturbative 
effects phenomenologically, which lies in the range: < k < 1- The equation of 
motion for the conformal mode is then obtained as 

- ^Bod^4> + Mle^^ {(ydl4> + 69„</)a^(/)} = 0. (3.3) 
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This equation comes from the trace part of the stress tensor. Here, there is no 
contribution from the Weyl and the matter terms, which are trivially traceless. 
The energy conservation equation from the time-time component of the stress 
tensor is given by 

{293</.a^^ - - 3M|e2^a^0a^(/. + e'^^p = 0. (3.4) 

Since the matter sector is conformal invariant, we here assume the stress tensor 
to be a form of the so-called perfect fluid with density p, although its contents 
may be different from that of ordinary matters. The contribution from the Weyl 
action also vanishes in the energy conservation equation. 

It is known that these equations have a stable inflationary solution for the 
region where the coupling constant tr is small. As mentioned in section 2, 
however, the coupling constant can not stay small in the inflationary space- 
time. We must take into account the gradual increase of the coupling constant 
along with the expansion. The inclusion of the dynamical effects is not simple, 
and the equations of motion become rather complicated containing non-linear 
as well as non-local terms. Bravely, we simplify it under the following physical 
consideration: the running coupling constant ir(p) is an operator which acts 
on the field (j), and it may fluctuate very much even on a smooth background. 
When the scale of local fluctuation of the order 1/Mp is small compared to the 
size of the system, 1/Aqg, we approximate the running coupling operator by 
its average under the spirit of the mean-field approximation. We rewrite the 
running coupling operator by time-dependent average replacing the physical 
momentum by 1 /r in equation (|2.8p : 

= /3olog(l/r2A|^)' (3-^) 

where r is proper time defined by dr — a{rj)drj with a — e'^. It shows that the 
running coupling diverges at the dynamical time scale, r = ta with ta — I/Aqq. 

Replacing the constant t'^ in the dynamical coefRcient (|3.2p by the time- 
dependent running coupling constant (|3.5p . and introducing the variable H = 
d{T)/a{T), where the dot denotes the derivative with respect to the proper time 
T, the trace equations of motion (|3.3p is written as 

^Bo(r) (^H +7HH + AH^ + ISiJ^i? + 6i/'*) - 3Af| [h + 2H'^^ ^ (3.6) 
and the conservation equation (|3.4p is 

-^Bo{t) (2HH -H^ + 6H^H + 3H^^ - 3A4H^ + p ^ 0. (3.7) 

For the region where the running coupling is small, these equations have 
a stable solution in which the scale factor exponentially grows up with the 
expansion time constant H ~ Hu, such that a(r) ~ 6^°"^, where 




HYy = \—Ah. (3.8) 
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As the coupling constant gradually increases, the value of the function Bg re- 
duces, and the value of H and its time derivatives increase. At the transition 
point T\ for < k < 1, the the third derivative of H diverges. In the case of 
K = I, the second derivative of H also diverges. The combination BqH, how- 
ever, vanishes in any case at the transition point so that the energy density is 
kept finite. 

The energy conservation equation implies that the density p starts almost 
vanishing in the inflationary era, and at the transition point it sharply increases 
to the finite value p(ta) = 3MpH^{T\)^ This behavior becomes more clear 
when we consider the time derivative of matter density, 

p + AHp = -^Bo(t) ( 2HH -H^ + 6H^H + 3H^) . (3.9) 

Initially Bq ~ 0, and the matter density is not increasing, while it increases 
sharply about the transition point where Bq changes drastically. 

Solving equations of the motion numerically, we obtain the solutions for 
H and p shown in figure 13. 1[ and a and p in figure 13.21 The simulation is 
performed starting at the Planck time defined by rp — 1/Hd and ending at 
TA = 1/Aqg. We choose three values bi — 7, 10, 15 for the coefficient of the 
Riegert action determined from the conformal field, where the values correspond 
to the Standard model and various GUT modclso The other parameters are 
rather indistinct because they depend on the non-perturbative dynamics of the 
traceless mode, and they are likely to be determined phenomenologically. Here, 
they are chosen as /3o = 0.6, ai =0.1 and k = 0.5. The ratio of the mass scales 
is set as Hd/Aqg = 60 and by this choice the number of e- foldings from the 
Planck time to the dynamical time becomes 

AA, = log4^ = 65. (3.10) 

This e-foldings gives the desirable number large enough to explain the evolu- 
tional scenario of the universe from the Planck time to the present and the 
primordial spectrum discussed in section 5. 

This simulation suggests that the quantum fluctuation of the conformal field 
covering over the whole quantum universe makes a transition to matter fluc- 
tuations at TA, and the expansion slows down from the inflating conformal 
space-time to the Einstein space-time as discussed below. The entropy is then 
generated to develop the thermal universe, which is regarded as the big bang in 
our scenario. The primordial spectrum we extract from the CMB observation 
is expected to be the reflection of quantum fluctuation of the gravitational field 
developed right before the transition. 

*This value is consistent with the radiation density for the Hawking temperature of deSitter 
space-time propotional to i?D , although the meaning of the temperature in quantum space- 
time seems obscure. 

5 For example, bi = 7.0 (A^a = 12, A^w = 45) for the Standard model, bi = 9.1 for 
SU{5) [Na = 24,Afw = 45), bi = 12.0 for 50(10) (A^a = 45), and bi = 17.7 for Eg {Na = 
78). 
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Figure 3.1: The time evolution of H and p for bi = 7, 10, and 15. iJo is nor- 
malized to be unity, and thus ta = 60 and inflationary solution for H becomes 
independent of bi. The dashed line denotes the Friedmann solution for H. 

Einstein phase Below the energy scale Aqq, the Einstein action becomes 
dominated, and the space-time makes transition to the Friedmann universe. 
Here, in analogy to the chiral perturbation theory as the low energy effective 
theory of QCD |33], we derive the low energy effective theory of gravity valid 
below the energy scale Aqg, and connect the solution in the inflationary phase 
with that in the Einstein phase smoothly at the transition point. 

In the case of QCD, dynamics of gauge fields disappears below the dynamical 
QCD scale, and meson and baryon become dynamical fields. In quantum gravity, 
although dynamics of conformal gravity disappear below Aqg, the metric tensor 
still remains as the dynamical variable in the Einstein gravity. As a consequence, 
the low energy effective action is given by an expansion in derivatives of the 
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Figure 3.2: The time evolution of a and p for bi ~ 7, 10, and 15 in a logarithmic 
plot, normalized by their values at ta {— 60) denoted by ca and p\, respectively. 
The dashed line denotes the Friedmann solution for p. 

metric field as 

ho^ = j d^x^{C2 + Ci + ---}, (3.11) 

where the lowest derivative cosmological term is not included because it is neg- 
ligible at the energy scale we are discussing. The two derivative term is the 
Einstein action 

C, = ^R + Cf^ (3.12) 

where is the conformally invariant matter lagrangian. 

In practice we employ terms up to the fourth order. Following the prescrip- 
tion of the low energy effective theory, £2 is used for calculating both tree and 
one loop diagrams, while £4 is used for the tree diagrams. The Planck mass is 
then regarded as the counterpart of the pion decay constant in chiral pertur- 
bation theory. The higher derivative terms are expanded by the inverse of the 
Planck constant, which can be justified owing to the ordering Mp ^ Aqq. The 
possible terms for £4 are given by 

p2 p2 p2 r> rpliv _}_rpilVrpM /n 1 o\ 

^ 1 -fi-^ii/j -"-/^i/Actj 7i.r2 -^Mi^^M I A/f4 \O.LO) 
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where is the stress tensor for conformally invariant matters. 

Since we restrict our effective action up to the fourth order, the Einstein 
equations of motion MpR^„ = is used to reduce the terms in £4. The 
equation of motion also imphes R = 0. Taking into account these equations 
and Euler relation to remove the square of the Riemann-Christoffel curvature 
tensor, the independent terms in £4 is simply given by 

Ci = j^R^^'R^,. (3.14) 

Here, the constant a is a positive value which is phenomenologically determined. 

The constant a acquires loop corrections coming from £2- The renormaliza- 
tion procedure is done by background field method [3 [S] about the background 
satisfying the Einstein equations of motion with the cutoff A (< Aqq). The 
constant a receives a finite renormalization dependent on the cutoff scale. The 
differences between two cutoff scales, A and Aqq, will be given by 

a(A) = a(AQG) + Clog(AVA2jG). (3.15) 

Here, ^ is a gauge invariant constant, because the Ricci tensor is divergence- 
free by equation of motion: V^R^v = V^T^ = 0. For the contributions from 
diagrams with internal matter loops, C, = {Nx + 3iVw + 12A'a)/120. Equation 
(|3.15p shows that the coupling a becomes small and the four-derivative term 
becomes irrelevant at low energies. It is worth commenting that this theory is 
valid at the low energy below Aqq, and therefore a ghost pole located at the 
order of Mp does not emerge. Therefore, the higher derivative action in the low 
energy effective action does not conflict with unitarity. 

We firstly solve the equation of motion for the homogeneous component in 
the low energy effective theory. It is given by the trace of the stress tensor, or 
the variation of the conformal mode, 

A4 (^H + 2H^^ +^{h +7HH + 4ij2 + UH^H^ = (3.16) 

for H as a function of the proper time. The energy conservation equation is 
obtained as 

- u4h'^ + P+£^ [-QHH + 3i/2 - ISH^H^j = 0. (3.17) 

When we restrict terms in the effective action up to the fourth order in 
derivatives, it would be valid at the energy scale sufhciently below Aqq- In order 
to fill the gap between Aqg and the low energy scale A, we naively interpolate 
the equation up to Aqg by assuming the coupling to be time-dependent: 

a(T) =ao-f Clog (^^) ^ . (3.18) 

where ao — a (Aqg). This form is given by replacing the cutoff with the inverse 
of the proper time in (j3.15|) . By the running coupling we expect to take into 
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account the space expansion effect as we have done for the coupHng tr{p) — > 
tr^l/r) in the last section. By this couphng, it is expected to vanish at low 
energy, and it is written in an non-negative form. 

Equations for both sides of the transition time ta are the same order, and 
we can smoothly connect the Einstein space-time with the inflationary space- 
time. The initial values of H, H and p are chosen to connect with those in the 
inflationary phase. The initial value of H in order to solve equation (|3.16p is 
given by solving the energy conservation equation (|3.17p . We can show that 
the Friedmann universe, H -\- 2H^ = and 3MpH^ = p, is a solution of equa- 
tions of motion (|3.16p and p.l7p , and general solutions are approaching to the 
Friedmann universe oscilating about this solution. The numerical simulation is 
carried out by setting the parameters as ao = 1 and C = 1- The results are 
given in figures [3T] and [3T2l in which Hd is normalized to be unity such that the 
solutions depend on the value bi. For a few steps in the period after the tran- 
sition, the space-time expands acceleratingly and the matter density decreases 
sharply during that period. It soon reduces to the Friedmann universe. 

4 Linear Perturbation Theory 

Since the inflation is a stable solution, it is expected that the amplitude of 
fluctuation of scalar curvature, SR, damps during the inflation, but does not 
vansish. The amplitude is estimated as follows: since scalar curvature has two 
derivatives, the size of the fluctuation at the Planck scale is estimated to be 
SR ~ TJq such that 6R/R ~ Hl,/12H^ ~ 10~^ in the inflationary background, 
where the denominator is normalized by the deSitter curvature with H = Hd. 
Near the transition point, the running coupling gets large and the energy scale is 
the order of Aqq. Then, the the scalar curvature contrast at the transition point 
is estimated to be 5R/R ~ K^q/12H'^. If we take the dynamical mass smaller 
than the Planck mass by two digits, the amplitude becomes the observed order 
of 10~^. This allows linear perturbations [Ml ISZj about the inflationary 

solution applicable from the Planck scale to the dynamical transition scale. 

The asymptotic freedom implies that the fluctuations of vector and tensor 
modes are relatively small compared to the scalar mode. In addition, the fluc- 
tuation we will discuss expands rapidly enough during the inflation to the size 
far from the horizon scale, and thus it is not affected by the dynamics near 
the transition point. The linear perturbation can be used for such a size of the 
vector and the tensor as well. 

Perturbations The conformal mode and matter density are perturbed about 
the homogeneous solution as 

0(77, x) 0(77) 4^1^(77, x), 

piv.^) = pin) + Spiv, (4-1) 

where (pirj) and ^(77) are solutions of the equations of motion (|3.3[) and (|3.4[) . In 
below, except in the appendix A, (j) denotes the homogeneous solution. 
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The stress tensor for the matter sector is traceless, which is taken to be 

rj^MA^ = 0, 



= e4'^|i(p + ^p + 4p^)<5,,+n,,j, (4.2) 

where Vi = diV + vj are velocity perturbations, where vj is the transverse 
component. The anisotropic stress is traceless, 11*^ = 0. The appearance of 
the conformal factor is due to the convention of the stress tensor defined in the 
appendix A. 

Hence, we consider the linear perturbations of field variables, (p{r], x), h^,^{r], x), 
6p{r],y:), and Vi{r],x). 

Gauge invariance Under the general coordinate transformations, S^gf^i, = 
QnX^uS,^ + guxV field variables are transformed within a linear approxima- 
tion as 

S^h^u = dij,^^ + d,y^^ - ^r]i^^dx^^ , (4.3) 
where ^fj, = r^nv^" . The traceless tensor mode is further decomposed as 

/loo = h, 

hoi = hj + dih', 

hij = h]^ + d^,hl' + ls,jh+[^-ls,,>jh", (4.4) 

where i,j = 1,2,3. hJ and hJ' are the transverse vectors and hj^^ is the 
transverse-traceless tensor. = d^di is the spacial comoving Laplacian. If 
is decomposed as ^° and — Sj + di^^, the general coordinate transformations 
are described in terms of ^'s as 

khj = dr,^J, 
5^hJ' = 2^J, 

khj;" = 0. (4.5) 
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For the matter sector, the perturbed variables are transformed as 



6^v 






= -d,C 








= 0. 



(4.6) 

The gauge invariant gravitational potentials so-called Bardeen potentials are 

defined by 

$ = if+^h-^h" + adn4>, 

* = - + ad^cj) + d^a, (4.7) 

where 

If we take the gauge h' ~ h" = 0, the Bardeen potentials are expressed as 
^ = ip + h/Q and = (p — h/2 such that the metric has the foUowing form: 

ds"^ = [- (1 + 2*) V + (1 + 2$) dx2] (4.9) 

for the scalar perturbations. The gauge invariant vector and tensor perturba- 
tions are defined by 

ri = hj-^dnhj' and /i^t. (4.10) 
For the matter sector, gauge invariant perturbations are defined by 

p p 
^ = ^+2^' 

Vi = Vj + ^dr,hj', 

= vj + hj. (4.11) 

Here, the vector variables satisfy the relation Ti + Vi — fli = 0. 

The gauge invariancc restricts the form of the stress tensor with the dy- 
namical factor for the conformal-field sector. In order to obtain the gauge 
invariant equations of motion, it is necessary to rewrite the function Bq by the 
field-dependent function B which transforms as a scalar under the coordinate 
transformation: 

5^B = eOxB = eOr^B. (4.12) 
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Using the a field, which is transformed as 5^a = the scalar function B is 
written as 

B^Bo-ad^,Bo. (4.13) 

The stress tensors for the Riegert action are given in the appendix A. Replacing 
the coefficient bi in the expressions by the dynamical factor biB, we obtain the 
equations of motion for the conformal-field sector. 

Linear scalar equations Let us write down the linear equations of motion for 
the gauge invariant scalar perturbations, in which the anisotropic stress tensor, 
Hij, does not contribute. We find four independent dynamical equations of 
motion for four independent scalar perturbations, $, 4', D, and V. We first 
consider the following two equations: 

T\ - 0, (4.14) 
^(TV-3^T,,j = 0. (4.15) 

These equations are independent of the matter sector, and we can obtain the 
Bardeen potentials by solving these equations. In terms of the gauge invariant 
variables, the trace equation (|4.14p can be written as 

+ (125^0 + 12a„0a^0 ~2f)-^} = 0. (4.16) 

Here, the field and the factor d^jBo in B are removed by using the ho- 
mogeneous equations of motion. The formulas in the appendix B are useful 
when equations of motion are transformed in terms of the proper time. The 
space-space component equation ()4.15p gives 

h ( A / 28 8 8 \ 

-|-M|e2'^{-2<I> - 2*} = 0. (4.17) 
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Equation (|4.17p is of a second order with respect to the time derivative. This 
equation plays an important role in connecting between the inflation and the 
Einstein phases. In the limit tr — s- 0, where the conformal field dominates, 
$ = is realized due to the vanishing of the Weyl tensor, while at the transition 
point where the coupling diverges, the configulation with $ = — 4' should be 
realized. 

The matter density fluctuation D is determined by the equation for (00)- 
component of stress tensor, and the velocity perturbation V are determined by 
the equation for (Oz)-component. We here consider the following two combina- 
tions: 

Too + 39,,(/.|^T,o = 0, (4.18) 
|^T,o = 0. (4.19) 

In terms of the gauge invariant variables, the energy conservation equation (|4.18p 
is written as 

Bo{t) ( (-2920 + 29,05,0 - I A + (25,30 - 45205,0) 5,$ 



87r2 



+d„cf, (-25^0 + 25,05,0 - 2f) 5,$ + (^-|^5,05,0 + ^ f<S> 

+5,0 (^25,20 - 25,05,0 + ^ 5,* + (-25^05,0 + 45^05^0) * 
2. ^ 2.2 



25,^0+ -5,05,0 
A^iq, _ 45,0^25,$ + + 45,0^2 5,* 



' <2(r) i 3^ ""^ " ' 3 
+M^e^'''2^^^+ e'^'f'pD = 0. (4.20) 

and equation (|4.19p is 

^Bo(r)|-^5,3$ -t- (^-^5,^0 ^5,05,0 5,$ - ^d,cf>f'^ 

+19,05^* + (^25,> - ^5,05,0 + 5,* + (^25^0 - ^5,0^?^^ vj/j 



tl{T) i 3^ " 3 
+Ml e^"* {25,$ - 25,0*} - I e'^'^pV = 0. (4.21) 

These equations involve up to the third order in the time derivative. Thus, 
substituting the Bardeen potentials obtained from (|4.16p and (|4.17|) into these 
equations, we can get the density perturbation D and velocity perturbations 
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V. Since terms from the Riegert and the Weyl sectors with dynamical factors 
disappear at the transition point ta according to the vanishing of the factors, 
and Bq, the values of D(t\) and V(t\) are determined from the value of 
^(■''a) (= — ^'('''a))- Thus, it is enough to compute the Bardeen potentials to 
achieve information necessary at the transition time ta- 



Linear vector equations The vector fluctuation is determined by the 
equation for (ij)-component of stress tensor, and the velocity perturbation 57^ 
is determined by the equation for (Oi)-component of stress tensor. We consider 
the following combinations: 

St., = 0, (4.22) 



To. = 0. (4.23) 



Equation (j4.22|) is written as 



hi 

+Mle^^\^dnT, + a„0T,| - 0. (4.24) 
and equation (|4.23p is 

+^M^e'''^fT,-^e^'^pn, = 0. (4.25) 

The second equation is composed of variables with at most the second order 
of the time derivative, while the first equation is of the third order. Thus, the 
value of rii(TA) is determined from the value of Ti(rA) as discussed before. 



Linear tensor equations The linear equation for the gauge invariant tensor 
perturbation is derived from the space-space component of the stress tensor, 
Tjj ^ 0, as 
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Bardeen Potential <I)(i[=10, m=0.05) 




Figure 5.1: The evolution of Bardeen potential $ for bi = 10 and m = 0.05 
Mpc"^ 

5 Primordial Power Spectra 

The observed CMB anisotropies are understood as a reflection of the Bardeen 
potentials on the last scattering surface [351 [3Z| ■ Before the universe was 
neutralized, we can trace the history up to big bang using the Einstein theory. 
The dynamics we described in the previous section enable us to trace back 
further to the Planck scale where the inflation ignites. 

The low niultipole components of the observed CMB anisotropies correspond 
to the fluctuation of the size of the order of the Hubble distance I/Hq, where 
Hq is the present Hubble constant. If the universe expands more than about the 
order of 10^*^ from the begining during the inflation and the Einstein periods, 
the fluctuation of the order of the Hubble distance originates from a fluctuation 
less than the order of the Planck length. The idea of inflation suggests that this 
order of the expansion is able to solve the flatness and horizon problems. Thus, 
we hope to observe the Planck scale physics from the CMB anisotropics. 

We solve the coupled system of equations, ((iTB)) . (|iT7)) . (g^ll), and (g^ll), 
numerically to obtain the values of the Bardeen potentials, the vector and the 
tensor perturbations. We start the simulation at the Planck time rp {— I/H-q), 
where the linear perturbation theory becomes applicable. At this time, the run- 
ning coupling constant is still small enough so that fluctuations of the conformal 
field are dominated by the Bardeen potentials satisfying the relation <t> = 5", 
and the vector and the tensor perturbations are relatively small. 

The scalar power spectrum is given by the two-point quantum correlation of 
Bardeen potentials computed by conformal field theory. We write the Fourier 
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Vector Perturbation (i);=10, m=0.05) 



Tensor Perturbation {/7/=10, «i=0.05) 




Figure 5.2: The evolution of vector and tensor perturbations, and hjj , for 
61 = 10 and TO = 0.05 Mpc^i. 

transform of a field / in a comoving momentum space as 

/(r,x)= J d3k/(r,k)fc-3g»k.x^ (5^^) 

and its two-point function as (/(r, k)/(T, k')) = \f{T,k)\^k^S^{k - k'). The 
initial profile of the scalar spectrum is assumed to be [39l [6] 

/ I \ "s — 1 

PpS(fc)^|$(rp,fc)|2=As(^-j . (5.2) 

Here and in the following we write k = |k|, and As is a dimensionless constant. 
TO = a{Tp)HD is the comoving Planck scale at the Planck time. The spectral 
index Ug is given by the anomalous dimensions of the scalar curvature, 



= 5 ~ 8 V / I ^ 1 ^ 2/61 + 4/6? + 0(1/6?), (5.3) 

1 - v/1 -4/61 

where 61 is given by (|2.6p . For the large value of 61, the index approaches to 
that of the Harrison-Zel'dovich 40 . 

The initial profiles of the vector and the tensor spectra are given by 

PpX(fc) = |T(rp,fc)p = Av(^^) 
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Figure 5.3: The scalar and tensor spectra at the transition point for bi = 7, 10, 
15 with m — 0.05 and m = 0.06. The initial amplitudes of bi = 10 for each m 
are set as y/Ag = 10~^ and \/Ar — 10~^. The others are arbitrarily chosen to 
be equal to the result of 6i = 10 at fc = m for each m. 




10"' 10" 10"' 10"^ 10"' 10"' 



t [Mpc"'] * [Mpc"'] 



Figure 5.4: The logarithmic plot of the scalar and tensor spectra at t\. 



where Ay and At are small dimensionless constants. The indices are given by 

Uy ^ nt ^ 0, (5.5) 

because the vector and tensor modes are conformally free fields with logarithmic 
correlations for the distance in the space coordinate. 

The spectrum implies that the fluctuation of the distance 1 /m at present 
originates from that of the Planck scale, I/H-d, at the Planck time. Beyond the 
Planck scale, the non-linear effects of the conformal mode will become signifi- 
cant. In order to include non-linear effects, we have to solve the model exactly 
or to use a non-perturbative method such as dynamical triangulation f41' for 
the regime beyond the Planck scale. 

The initial values of the Bardeen potentials are set as $(tp, k) — ^'(tp, k) = 

^Pp[(fc) at the Planck time. Since the Bardeen potentials are fluctuations about 
the homogeneous solution, the boundary values at A: = satisfy the Dirichlet 
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conditions 

$(r,0) = ^'(t,0) = 0. (5.6) 

Since the running coupling diverges at ta, equation (|4.17|) implies that the 
Bardeen potentials at the transition point satisfy the relation 

$(rA,A:) + ^'(TA,fc)-0. (5.7) 

We numerically solve the coupled equations (|4.16p and (|4.17p imposing the 
boundary condtions, (|5.6p and (|5.7[) . The initial conditions for the time deriva- 
tives of $ up to the third order set to zero and satisfies the boundary (|5.7p 
at Ta. The initial values of the vector and the tensor perturbations are set as 
T,i(Tp, k) = \/Ay and hJ^{TY>, k) — \J and the time-derivatives of these fields 
are also chosen to vanish at rp. 

The density and the velocity perturbations are computed using the equations 
(|4.20p and (|4.21[) , respectively. The values at the transition point are calculated 
using the Bardeen potential to be 

oH(ta) m'^ 
kVU,k) = l^!R^e-^^-U{T^,k) + H{TAMTA,k)), (5.8) 

where J\fe is the number of e-foldings defined by (|3.10p . Here, we multiply k on 
V in order to make the expression dimensionless. Also the vector perturbation 
fli is computed using equation (j4.25p . resulting 

To derive these values, we use the value p{t\) = 3MpH{ta)^ ■ Due to the 
exponential factor, these values are small enough to solve the flatness problem. 

Let us examine the parameters to be used for the simulation. The results 
of simulation are given in figures 15.11 and 15.21 The value of the coefficient in 
front of the Wess-Zumino action is taken to be bi = 10, which gives ng = 1.25, 
and the values of other parameters to determine the homogeneous background 
are the same as those used in section 3. The comoving Planck scale is taken as 
m = 0.05 Mpc~^. Since Hd is given by the order of 10^^ GeV, the scale factor 
at the Planck time a{Tp) becomes the order of 10~^^, where the present scale 
factor is normalized to be unity. On the other hand, after the transition occured 
at the energy scale Aqq 1.1 x 10^^ GeV, where Mp = 2.4 x 10^^ GeV, the 
universe expands about 10^^. This implies that the space-time expands about 
the order of 10^° in the inflationary period starting from the Planck time and a 
few accelerating expansion period right after the transition before reducing to 
the Friedmann universe. The number of e-folding for these periods corresponds 
to 70. It is consitent with the solution given in section 3. The initial amplitudes 
of the Bardeen potentials, the vector and the tensor perturbations are taken to 
be y/AE = 10-^ and ^/AU = VA^ = 10"^ 
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The amplitudes of the Bardeen potentials reduces much slower than the 
damping rate during the inflation. Existence of the gradually varing component 
can be seen in the analytical solution of approximate equation of motion valid 
in the intermediate stage of inflation as shown in the appendix C. The tensor 
fluctuation is preserved to be small, while the vector fluctuation increases near 
the transition point. However, in the Einstein phase, the scalar and tensor fluc- 
tuations survive until present epoch, while the vector fluctuation soon disappear 
and we cannot observe it. The primordial scalar and tensor spectra evaluated 
at the transition point are deflned by 



The results are given in figures [5751 and [5T^ where we plot the cases of bi{ns) = 
7(1.41), 10(1.25), 15(1.15) with m = 0.05 and 0.06. The simulation is carried 
out for the region k < 0.15, using the Fortran software, BVP -SOLVER |42j. 

It seems that the linear approximation is not applicable for the high momen- 
tum region k"^ /ni? ^ 1, because such a fluctuation corresponds initially too far 
beyond the Planck scale so that there appear coefficients with large values in 
the equations which violate the applicability of the linear approximation. The 
patterns of the spectra are sensitive to the values of m and 6i, while the pa- 
rameters determing the dynamical factor Bq are insensitive to that, apart from 
the magnitude of amplitude. The patterns of spectra are also preserved in the 
Einstein era, although the amplitudes may slightly change in a few moment 
after the transition. 

The initial spectra we employed here do not cut low momentum components 
required to explain the observed suppression of low 1-components in the angular 
power spectra. In order to explain it, we meight need to seek a reason in 
the ambiguity of the choice of initial spectra which should reflect original non- 
perturbative dynamics of the traceless mode for the long-distance two-point 
correlation functions. 

6 Conclusion 

We have constructed an evolutional model of the universe based on the con- 
formal gravity with the symmetry breaking Einstein action. The universe starts 
from a quantum state at the Planck time and grows up exponentially. As a con- 
sequence of the asymptotically free dynamics, there appears a strong coupling 
phase where the structure of the space-time changes drastically. The transition 
takes place at the dynamical energy scale of gravity, where quantum space-time 
with conformal invariance makes transition to classical Einstein universe. Then 
the universe grows up to the present size, such that the size of the Planck length 
at the Planck time extends to the size of the order of 10 Mpc distance today. 
We have suggested that at the transition, field fluctuations freeze to localized 
objects, and they eventually decay into the classical matter driving the universe 
into the big bang phase. 



pS(fc) 
P\k) 



l<i>(TA,fc)P, 



(5.10) 
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It was shown that the linear perturbation is applicable on the inflationary 
background for the momentum range which covers the size of fluctuation ob- 
served as CMB anisotropics today by COBE and WMAP. The evolutions of 
gravitational scalar, vector and tensor fluctuations have been evaluated. Since 
the initial fluctuations are provided by conformal fleld theory, we expect that the 
amplitudes of tensor and vector fluctuations are relatively small in comparison 
to that of scalar fluctuations. The scalar fluctuations are getting small during 
the inflation, and the tensor fluctuation is preserved to be small, while the vector 
fluctuation is getting large near the transition point. However, since the vector 
fluctuation disappears in the Einstein era, the tensor fluctuation in addition to 
the scalar fluctuation may contribute to the primordial spectra constituting the 
observed CMB anisotropics. 

The condition Mp ^ Aqq is significant to make the inflationary scenario. 
It also implies that quantum effects turn on the size much larger than the 
Planck length so that not only the space-time singularity but also the horizon 
of an elementary excitation with the Planck mass disappear. Furthermore, the 
strong repulsive effect in quantum gravity which causes the inflation also erase 
the singularity inside a black hole by balancing the pressure of the collapsing 
matter. However, if the collapsing goes too far and exceedes the balance, the 
black hole may explode something like as mini-inflation. 



A Stress Tensor for Each Sector 

In the appendix A, the conformal field (j) denotes a full field including the 
homogeneous and non- homogeneous parts, such that = 0(77, x) (|4.1[) . To 
obtain the linear perturbation theory, the stress tensor has to be expanded by 
the fluctuating field ip. 

Equations of motion are obtained by the variations of the effective action as 



where the indices of stress tensor T^^ are contracted by the flat background 
metric rj^,^. The ordinary stress tensor defined by the variation of the physical 
metric gp,i, is denoted by T^,^. Furthermore, we write the stress tensor defined 
by the variation of the metric g^^ as T^^. Indices of these stress tensors are 
contracted by g^i, and g^i,, respectively. The difference of T^^ and T^i, is just 



Appendix 



6T 




(A.l) 
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a conformal factor such that T^^ = g-^'^T'^^ and T^'^ =_ e'^'^T'^, and so on. 
Within hnear perturbations in h^^^, the relation between T^^ and T^^ is given 
by the symmetrized product T^i^ = 'rjx[^T^^y 

The density and velocity perturbations in the matter sector are defined by 
the physical stress tensor on the metric g^n as 

T^l = -{p + 5p), 
T^o = -{P + P)v\ 
T^" = {p + P){v,+h,,), 



T^) = {P + 6P)S'^+W^. (A.2) 



where Vi — (jijV^ . For conformal invariant matters, the equation of state is 
given hy w = P/ p = 1/3, and SP/Sp = 1/3. The trace of the anisotropic stress 
vanishes, Il'j = 0. The stress tensor is transformed as S^T^i^l — d^S^^T^^-^ ~ 
dx^'^T^^ + S/^dxT^'^l under the general coordinate transformation. 
Hence, the dynamical equations of motion are given by 

T,, = T«„ + Tjr. + T™ + = 0, (A.3) 

where R, W and EH denote the stress tensor for the Riegert action, the Weyl 
action and the Einstein action, respectively. They are given below, while T^^, 
is given by (|4.2p . 

Riegert action 



T 
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+Ah''^dxd^d^,,(l)d,)^ + Ad(^,,h^^dxd^(f>d,)^ + Ax^'dxd^^^d.jcj) 

+lh^^d^dx<f>d,d„(P + ^d^^h^''d,)dx^d^^ + Ad^h^d.^dx^da^ 

+ld''h^,dxd^<l>d''cf> + U^d,h^''dx<l>d,cl> - Ad^d(^hl^dxcl>d„^ 



-2d^d''h^,dx^d„4> - 49^x(^9,)(/.aA0 + Auhl^d,)4>dx4> 

'^d(^Xu)dx<pd^(t}~- '^^h^ydx(t)d^(t)+ ^dxx^d^j,<l)dy<j) 
-Ahl^d,)dx4>0^ + 2h\^d,)U^dxct> + 2/iA(^9^n</)a.)</> 

-^h\^d,)d''cj,dxd,(l> - ^h\^d,^dxd,<i)d''ci, 
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2 



^h^'^dxd^d^d,^ 



^d^d,h^^dxd„<l) - ^x^dxd^d,^ 



-d(j.h^''d,)dxd. 



+Ad%^hZ^dxda 



2d^d''hnudxda 
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3 d(^u,Xu)^(t> + ^Oh^^^Oc/) 



-2dxx^d^d,4> + 45aX(m5.)5V - ^ahld,-)dx<^ + -dxd^^x^d. 



.-hld,)dxUcl> 



-Id^^dxdack + l^X^dx<l> + 2dxx''^<^ - Id^d^x" dxct> 



■\d^,d,dxx^ + \vt^u^dxx^ 



(A.4) 



where Xju = d\h\ and □ = d^dx = —9^ + The trace of this tensor is given 
by 



= ^|-2n2<^ + 4/i'^''na^S,</. + 4d^h'"'dxd^d,<p + Ax^dxOcp 
+ADh^"'d^d,cP + 2Dx^dx<l> + ^dxX^O<i> - ^d^d^x'^dxcp 



+ lodxX^ 



It is also derived by the variation of the conformal mode. 
Weyl action 



(A.5) 



1 



iA fiCT ^linear 



= -| - 2n5(^x.) + \d^d,dxx^ + ^'?^.°9aX^| • (A.6) 



This tensor is traceless: T^'^ 



0. 
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Einstein action 

-5(mX.) + l^h^. - 2/if^S,)aA0 + 2hl^d^^(f>dx<t> 

Q^AX^ + 2h^-dxd,<f> + h^^dx<Pd,cj^ + 2x^dx<p^ I . (A. 
The trace of the stress tensor is given by 

= Mle^-^j-en^ - Qd^(pdx(j) + dxx^ + Qh^^Oxd^cj) 

+6x^5a</) + 6h^''dx<^da(p} . (A. 

B Proper Time Representation 

The proper time is defined by the equation dr = a{T)dr], where a(r) = e'^^ 
is a solution of the homogeneous equation of motion. Then, we obtain 

dri = adr , 

d'^ = a^df + Hdr), 

5^ = a^{d^+3Hdl+(^H + 2H^)dr}, 

= a^{d^ + 6Hd^+(iH+nH'^^d^+(H + 7HH + 6H^^dr} 

(B. 

and 

dn<p = aH, 

dl(j) = a'^(H + H^^ , 

d^(j) = (if + 4HH + 2H^^ , 

d^(f> = a'^(H+lHH + AH^ + lB>H^H + QH'^^, (B. 
where II{t) = a{T)/a{T). 
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C Analytical Study of Linear Scalar Equations 



In this appendix, we study solutions of the coupled linear scalar equations, 
(|4.16p and (|4.17p . analytically. We here simplify the equations as follows: the 
coupling tr is a small constant such that the Hubble parameter are constant 
normahzed to he H = Hd/^/Bq = 1 and T = biBotr/Sn^ <C 1. Furthermore, 
we neglect the momentum dependence. Such a situation will be realized in 
several times after the inflation starts because the momentum dependence, k/a, 
becomes negligible as the scale factor grows. Then, we obtain the following 
equations for the scalar perturbations: 

-2 $ -14 $ -36$ - 48$ + 2 ^ +14* + 36* + 48* 

fi) + 4$ - * - 4*) =0, (C.l) 



+6 



4 •• 16 • 20 4 • 4 

-$ H $ H $ * + -* 

3 3 3 3 3 
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+ -2($ + *) = 0. (C.2) 

Introducing the variable / = * — $, these equations are written as 



/ +7/ +15/ +12/ = 0, (C.3) 

$-(i + ^T)<i>-lr<i> = -/-(i + ^t)/-1t/. (C.4) 

The first equation has the solution 

/ = ciC"''^ + C2e-5^sin ^^T^ +C3e"2^cos ^^rj . (C.5) 

Substituting this solution into the second equation, we obtain 

$ = (ai+ci)e-^ + (a2+C2) (^l-^Tr^ +(a3 

, 360 - 7r , V3c2 + 5c3 
+^^^80^' .^cos 

5C2 - VS cs 
14 





e--sin V^ . (C.6) 



Here, we expand the time constant in T and retain within the first order. 

Since we now consider the exactly zero-momentum case, this solution in- 
cludes both the homogeneous vacuum mode and the fluctuation mode we seek. 
At T = 0, the vacuum mode $ = * = w satisfies the following equation: 

'oj' +6 'uj +8ui - 3d; - 12a; = 0, (C.7) 
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which is obtained from equation (|C.1|) , while the second equation (|C.2p is trivial 
in this case. This equation has an inflationary mode, e"^, and three decaying 
modes, e^'^^ , e-^''/'^ sm{V3T/2) and e'^'^/^ cos(\/3t/2). Thus, by definition, 
the homogeneous vacuum modes, which behave like these at the T = limit, 
must be removed from solution (jC.6p . Thus, the fluctuation, <I>, behaves for 
T< 1 as 

$ _ 1 _ Ltt, (C.8) 
12 ' ^ ' 

where the exponential damping mode is neglected. This mode appears at the 

intermediate region in time evolution where the momentum dependence becomes 

negligible (see figure [5TT|) . 
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